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ABSTRACT 

A sizable fraction of the total energy density of the universe may be in heavy particles 
with a net dark U(l)' charge comparable to its mass. When the charges have the same 
sign the cancellation between their gravitational and gauge forces may lead to a mismatch 
between different measures of masses in the universe. Measuring galactic masses by orbits 
of normal matter, such as galaxy rotation curves or lensing, will give the total mass, while 
the flows of dark matter agglomerates may yield smaller values if the gauge repulsion is 
not accounted for. If distant galaxies which house light beacons like SNe la contain such 
dark particles, the observations of their cosmic recession may mistake the weaker forces for 
an extra 'antigravity', and infer an effective dark energy equation of state smaller than the 
real one. In some cases, including that of a cosmological constant, these effects can mimic 
w < —1. They can also lead to a local variation of galaxy-galaxy forces, yielding a larger 
'Hubble Flow' in those regions of space that could be taken for a dynamical dark energy, or 
superhorizon effects. 
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1 Introduction 



A variety of cosmological observations has revealed the prevalence of a dark sector in our 
universe. Typically we model it as a dark energy component, comprising some 73% of 
the total energy contents, and a dark matter component, contributing another 23% of the 
energy, with the rest in the form of familiar baryons (see, e.g., [1]). The simplest models 
of dark energy and dark matter posit that they are just two separate entities, most often 
a cosmological constant and a WIMP particle, in stark contrast to the visible sector which 
is known to contain a plethora of stable modes. This may at least make us curious that 
whatever birthed our universe may have been just as creative in populating the dark sector, 
as it has been in the visible sector. 

Many of our theories would appear to suggest this possibility. The additional dark 
modes are not yet necessitated by direct observational evidence 1 . However, promising micro- 
physical theories such as those which arise in string theory, for example, admit a dark 
population of the universe which may be at least as diverse as the visible one. There are 
many candidates for dark matter particles with a variety of interactions between them, 
with large masses and weak direct couplings to the visible matter. Among them some may 
indeed behave just like the WIMP, and be a dominant component of the relic dark matter 
abundance. In this work we shall consider what happens if it is not alone, focusing on 
signatures of additional dark matter with long range interactions, which for simplicity we 
model by a Maxwell-like theory. Various aspects of such theories have been studied in the 
literature [3]- [11], where dark matter was taken to carry either visible or dark charges, and 
where 'dark electrodynamics' have been dubbed U(l)' theories. They are not excluded as 
long as the dark charges are small enough to avoid too much cooling of virialized dark 
matter in galactic halos, which prevents their collapse into the galactic plane and preserves 
the observed halo structure. For TeV-mass particles, the bound on the dark fine structure 
constant is a' £ 10~ 3 , with stronger couplings for heavier particles, allowing for quite a large 
parameter space for such theories [11]. 

Here we are interested in extremely weak dark matter couplings. From the phenomenolog- 
ical point of view, clearly there is no limit as to how weak these couplings might be. However, 
theoretical arguments based on the considerations of UV completions of theories with weak 
gauge forces, namely holographic considerations and charged black hole (in) stability, point 
to the conclusion that in a given effective field theory regime, gauge forces should not be 
weaker than gravity [12]. These arguments also suggest that if the gauge theory coupling 
strength is fixed at g, its effective description should break at a low energy cutoff A < gMpi. 
Furthermore, there ought to be stable particles lighter than A, which obey the inequality 
M <, QMpi, and may saturate it if they arise from some fundamentally supersymmetric 
theory. 

In what follows, we will take the U(l)' couplings to be of gravitational strength, as is 
frequently the case in string theory. This yields g ~ 10~ 15 , and so by the arguments of [12], 
it implies the existence of a {7(l)'-charged lightest stable particle in the dark sector, with 
mass m < A ~ 10~ 15 M P i ~ TeV. If this particle also has short range weak scale interactions 

1 See, however, [2] for an interesting suggestion. 
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with the Standard Model and/or the remaining dark sector, by the arguments of [13], it may 
end up with the right relic abundance to contribute to the dark matter in our universe. 

If such particles inhabit our universe, they would have a number of interesting cosmologi- 
cal consequences. To leading order, some of them would be similar to those already explored 
in [14, 15, 16], where an additional Yukawa force due to a scalar exchange was assumed to 
operate in the dark sector. The ensuing bounds on the coupling strength, originally quoted 
to be of the order of gravitational coupling [14], were recently pushed down to about 10% 
of gravitational strength [16], with the underlying assumption that all of the dark matter 
experiences this extra force. Clearly, the bounds will be correspondingly weaker if only a 
fraction of the dark matter carries dark U(l)' charges. In this work we will explore what 
happens if all the U(l)' charges have the same sign 2 . We should note that the consequences 
of a weak unbroken U(l)' in a neutral universe have been examined recently in [11], where 
the authors considered interesting questions of dark chemistry and plasma dynamics 3 . In 
the case we consider, with feeble U(l)' of dominant strength, the effects one expects should 
be primarily competitive with gravitational phenomena. The cancellation between their 
gravitational and gauge forces may yield a mismatch between the masses measured by U(l)' 
singlets, such as baryon aggregates, and £/(l)'-charged objects, such as distant galaxies that 
may contain charged dark matter. Indeed, normal baryonic matter will follow orbits that are 
specified by the gravitational field alone, enabling us to correctly determine the total mass 
of the source, as is the case when studying galactic rotation curves or lensing. In contrast, 
the motion of charged dark matter agglomerates is controlled by both the gravitational field 
and the dark gauge repulsion. If one does not account explicitly for the gauge repulsion in 
the orbit determination, one would underestimate the mass of the charged source. 

This phenomenon directly affects the toils of an observer who extracts the Hubble rate 
from the cosmic recession of distant galaxies, where the standard candle SNe la reside along- 
side charged dark matter particles. If the observations are calibrated to the same value of H 
obtained from, e.g., CMB or other visible matter observations, our observer may mistake the 
weaker net forces for an extra 'antigravity'. This could lead to an interpretation involving ef- 
fective dark energy equations of state smaller than the real one, including mimicking w < — 1 
without any ghosts, phantoms or other similar pathologies. It could also cause a mismatch 
between the measurements of the dark matter density fraction CIdm using purely baryonic 
probes of dark matter and probes that track the evolution of total dark matter distributions, 
which may contain charged dark matter particles. The former are, of course, to leading or- 
der blind to dark matter charges, whereas the latter are sensitive to dark gauge repulsion. 
Further, galaxy-galaxy forces may change locally if the density of charged dark matter varies 
from place to place. This will affect recession velocities, causing deviations from the mean 
'Hubble value' that could be taken for a dynamical dark energy. Such dynamics might be a 
useful model to explain recent claims of the observed bulk flows in [17]. 

2 This would clearly require some early universe dynamics to explain a net non- vanishing charge in our 
borough of the universe, such as a 11(1)' charged inflaton condensate. 

3 It is, of course, rather difficult to imagine that such issues will be very important in our case when 
the binding force is as weak as gravity, since the binding energies in the 'dark sector' in our case will be 
g 4 ~ 10~ 60 times smaller than those in the normal sector, so that the bound states just won't be very bound. 
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2 Faradayan cosmology 



To demonstrate our claims, we start with a simple, yet sufficiently illustrative argument 
involving a derivation of the Hubble equation using Newtonian cosmology with charges. In 
other words, a Faradayan cosmology! This is in fact precisely the same approach that Lyt- 
tleton and Bondi used in 1959 to argue that a small fractional charge difference between an 
electron and a proton can yield Hubble expansion in the Steady State universe [18]. Indeed, 
if the charges of the electron and the proton had not cancelled each other exactly, hydrogen 
atoms would have had a small monopole charge q. Let us consider what happens if we take 
the cosmic matter to be charged under some U(l) field, with all charges having the same sign. 
Assuming for the moment that they fill the universe uniformly 4 , we model the extragalactic 
recession by the Hubble law v = HR, where R is the distance between the observer and a 
receding beacon. Describing the effects of gravitational and gauge forces by Gauss's theo- 
rem, such that the contribution to net acceleration by sources outside of the sphere centered 
at the observer's location vanish, the net contributions to acceleration comes from interior 
sources, and are given by Newton's and Coulomb's law, — G^ 22 ^ and respectively. We 
can now use total energy conservation to relate the Hubble parameter H to the matter 
contents inside the sphere, which is known to reproduce the correct General-Relativistic 
Friedmann equation. Taking the distant beacon to be marginally bound, with total energy 
E = corresponding to a spatially flat universe [18, 19, 20, 21], a little straightforward 
algebra yields 

2 8ttG n 8vr/gx /?n 

n — ^ Ptotai T\ , ~ , , , Pcharged , UJ 

O O ^Tfl' probe \ ff\ ' background 

where ptotai is the total cosmological mass density, p c harged is the mass density of charged 



particles, ( — ) is the charge-to-mass ratios of the distant beacon whose recession velocity 

\ m / probe 

is probing the cosmic contents, and ( — ) is the mean charge-to-mass ratio of the 

V m / background 

charged fraction of cosmic mass. 

Note that the correct interpretation of this equation is as follows: an observer would 
measure H, extracting it from the recession velocity of a distant object. If the object were 
neutral, one would simply get the true ptotai of the universe. If the object were charged, on 
the other hand, and the observer did not account for that effect explicitly, one would need to 
increase ptotai, in the regime of distances probed by the observation, relative to what would 
have been seen by a neutral probe in order to get the same value H . Said in yet another way, 
recall that the Hubble velocity measures the velocity of a particle that is marginally bound 
to a system of masses, at a location y from the observer who is pinned to a fixed location in 
the mass fluid. This means that if a probe particle is neutral it must have greater velocity 
in order to escape than a charged particle at the same distance. Essentially, the Hubble 
velocity is the analogue of the Newtonian escape velocity, and it will be smaller if some of 
the attracting masses are screened by the gauge charges. 



4 Here we ignore the paradoxes one faces in dealing with uniform distributions of charge in an infinite static 
space, with wisdom after the fact that such issues may be naturally resolved in time-dependent spacetimes 
such as an expanding universe, where we can treat a uniform distribution of 'charges' (a.k.a masses) in a 
perfectly meaningful way as an expanding universe. 
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Now, clearly, iff- 21 ) f- 2 -) ^Gn, and if the charged particles were a significant 

\m / p ro b e \ m J background 

fraction of the cosmic contents, the gauge repulsion would affect the energetics of the universe 
and influence the expansion rate. Had the typical matter in the universe been the hydrogen 
atoms, with masses m ~ GeV, the gauge repulsion would have significantly affected the 
cosmic balance for net monopole charges q ~ m/Mpi ~ 1CT 18 , as noted in the theory 
of Lyttleton and Bondi [18]. We can now recognize that this idea was really stemming 
from the mystery of nature which we have since identified as the gauge hierarchy problem: 
the tremendous disparity between the strengths of gravitational and electromagnetic forces. 
In that day, this at least motivated people to seek ways to improve the bounds on q in 
terrestrial experiments, culminating in the results of King [22], who set the bounds on the 
electromagnetic fractional charge difference between an electron and a proton to q £ 10~ 22 , 
rendering it cosmologically irrelevant and laying the Lyttleton-Bondi theory to rest for a 
time. 

A simple way to elude the King bounds and resurrect long range gauge theory cosmolog- 
ical effects is to take the dark sector to be charged under some U(l)'. Indeed, as discussed 
in [11], the dark sector U(l)' couplings can be significantly stronger than the King bounds: 
q' 2 ~ a' < 10~ 3 as opposed to <, 1CT 44 . As already stated in the introduction, we now 

take our charged U(l)' mode to be a TeV mass dark matter particle carrying integer units of 
dark charge q' ~ g ~ 1CT 15 , easily satisfying the dark sector bounds. If a significant fraction 
of dark matter is charged in this way, with similar fractions in extragalactic space and in 
galaxies and clusters, the gauge repulsion will affect the Hubble expansion as measured by 
the recession of faraway objects, because it will cancel the attraction between the charged 
masses. 

Indeed, consider Eq. (1) applied to the case where the universe is filled with a cosmolog- 
ical constant, neutral heavy dark matter and some charged heavy dark matter (with 'heavy' 
merely meaning, nonrelativistic particles now, with masses >, TeV). At some redshift 
z = 1/R(t) — 1, where R(t) is the cosmological scale factor (playing the role of the distance 
between an observer and a receding beacon) we have ptotai = A + p ne utrai + Pcharged, with 
Pneutrai and p c harged both scaling as inverse spatial volume: p(z) = p(0)(l + z) 3 . Substituting 
this in the formula for H given in Eq. (1), we find 

H\Z) = ^ A + ^{ Pncutral (0)+[l- ^H^j ^backgrounds) ] p charged (0) | (l + z) 3 . 

(2) 

This shows how the charged dark matter contributions get suppressed by gauge repulsion 



as a function of distance, or redshift. The explicit redshift dependence in (^J ^ (z) and 
( — ) (z) is not of microscopic nature, but simply notes that an observer utilizes 

\ m J background 

intrinsically different measures of cosmic expansion over the cosmic distance ladder, with 
different charged matter contents. This is the key point in what follows. For example, 
note that at very large redshifts one could extract the information about H(z) from CMB 
or baryon distribution variation, which to leading order are not directly affected by U(l)' 
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repulsion 5 . Since these probes are not charged under U(l)', we take (jL^J ^ (z) = 0. In 

contrast, at low redshifts one would use the luminosity- redshift diagram of Type la SNe, 
which reside in galaxies that may contain charged dark matter, and so feel gauge repulsion. 
In this case, our probe is charged, so ^ (z) will be non-zero. Even more prosaically, 

it may occur that as cosmic structure begins to form from adiabatic uncharged primordial 
seeds, the fraction of charged dark matter in the collapsed structures will initially grow, 
having started from zero, as long as the linearized approximation is valid and attraction is 
due to the primordial seed fields. The 'absorption' of charge will be quickly saturated by 
nonlinear effects, when the repulsion from the matter already in the structure diminishes the 
net force acting on the charged dark matter, while the mass will continue to increase, and so 
at very small redshifts one would expect a decrease of q/m. Clearly, this would also give rise 
to some non-trivial z dependence in the background charge to mass ratio, ( — ) (z). 

V rn / background 

The task of determining the details of the evolution of q/m is difficult, but nevertheless very 
important, and we plan to address it in detail in future work. 

If one normalizes Hq by large z observations and re-expresses (2) in terms of fractional 
dark sector densities fij = — , where p cr {z) = 3 ^q^ is the critical energy density of the 
universe, one finds 



(™) robe^ (™) back round ^ 
tt\(z) EI2 ^ ^ baCkgr ° Und fichargcd(^) + ^DM(^) • (3) 



Now suppose that using a set of fixed probes, which may experience extra forces, one does 
not account for the presence of dark charges and their repulsion, but merely identifies their 
contributions as 'not dark matter', lumping them together with the cosmo logical constant 
into 'dark energy' according to 



(z)( ± ) (z) 

n DE (z) = Slx(z) - N "* /pr ° bc ^ yback ^ un d fi chargcd (^) . (4) 

Applying this to the low z measurements shows that the gauge compensation of charged 
matter contributions effectively decreases the dark energy contribution Q\, and so to fit the 
data at lower redshift one may to increase the phenomenological value of A by hand at low 
z. That would change the phenomenological value of at high and low redshift. This can 
contaminate the extraction of the effective equation of state of dark energy from the data, 
since it may simulate the effect of increasing dark energy with decreasing redshift. Indeed, 
using the phenomenological definition of total dark energy f2 DE in Eq. (4), we can obtain 
the effective equation of state, wbe es, by using the definition of the equation of state of any 
given component 1 + Wi — —3^- = |(1 + z)d z ln(f2;iJ 2 ). So using Eq. (4), we find that the 



5 Of course, the neutral particles will be indirectly affected by the charges, by the fact that the structures 
in the universe, through which the neutral particles propagate, will be influenced by the extra forces. That 
does not affect to leading order the calibration of the luminosity-distance relationship, aside from the effect 
of inhomogencities which can be dealt with in the familiar manner, ft will however provide additional 
signatures of charges and forces, which should be considered separately. 
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effective dark energy equation of state is given by 



n (z)M (z) 

-, a 'charged^/ \ m J probe V m / background / r \ 

w »- = - 1 --iw3 <r„ ■ (5) 

where we have substituted w\ = —1, and u> cna rged = 0. In the derivation of this equation 
we have ignored ^-derivatives of ( — ) (z) ( — J (z) since that dependence does not 

\ m / probe V m / background 

really represent an intrinsic variation of the background, and merely reflects a change of the 
geometric probe at different redshift, as we have already explained. It is now evident that 
the misinterpretation of gauge repulsion as a dark energy effect could mimic u>de c ff < _ 1- 
This possibility has recently attracted considerable theoretical thought (see e.g. [23, 24]), 
and spurred much exploration of data (e.g. [25]). 

We have - for simplicity - assumed the 'real' dark energy to be a cosmo logical constant, 
as noted above. If the 'real' dark energy weren't a cosmological constant, but had equation 
of state, %e reah Eq. (5) would read 



O MM (z)W (z)\ 

I , . s 'charged \ z ) V W probe v ' V m J background v ' 
»DE off — %E real | 1 H 77 ~, \ ^ 



N 



(6) 



Note that this immediately follows from Eq. (4). Indeed, as the universe expands and 
the charged dark matter density redshifts away, the misplaced negative term in Eq. (4) 
becomes smaller, and so the 'net' dark energy Qde appears to grow even faster. We stress 
again that the key for this effect is clearly the use of different probes of cosmic expansion 
at different z, which have intrinsically different ( — ) , such that one may misinterpret 

the gauge repulsion as extra 'anti-gravity'. Without this one would merely renormalize the 
dark energy and dark matter contributions to the Friedmann equation relative to their mass 
fractions in the universe, as is clear from the equal scaling with z of both the neutral and 
charged dark matter contributions, for a fixed value of f- 2 -) f- 2 - ) . The point 

V m / probe V m / background 

is that gauge repulsion contributions scale the same way as neutral dark matter ones, as 
is manifest in Eq. (2). Without at least some change of the charge-to- mass ratio of the 
probe of the expansion and/or the background, the gauge repulsion alone could not affect 
the effective equation of state of dark sector 6 , or support a mirage of cosmic acceleration for 
long enough. 



3 From General Relativity to Faradayan cosmology 

The Faradayan cosmology description which we have used to write down Eq. (1) may appear 
dubious. After all, the standard cosmological model is firmly grounded on General Relativity. 

6 In their 1959 paper, Lyttlcton and Bondi wanted gauge repulsion to drive the growth of cosmic scale 
factor. To get it, they were assuming a constant charged particle number density, resorting to the framework 
of the Steady State universe, where matter is continuously produced out of nothing. Without this assumption, 
the redshift of number density precludes the change of expansion rate, as we have seen here. 
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The Newtonian approach however captures the essence of its physics correctly [19, 20, 21]. 
In the case of interest to us, this can readily be seen from an exact solution of Einstein's 
equations describing a system of charged massive particles in any FRW environment, where 
the charged particle masses are exactly compensated by their gauge repulsion. It is a gen- 
eralization of the McVittie solution [26] which includes any U(l) charges [27], and in the 
extremal limit where the charge is equal to the mass, it extends to a multicenter configura- 
tion thanks to a 'no-force' condition for stationary particles with q/m = ^/Gn [28]. Taking 
the cosmological spatial sections to be intrinsically flat, the solution is given by [27, 28] 

ds 2 = -d^ + R 2(t)n 2 (t,£)d&, 



n 2 (t,x) v ' v ' ' ' n(t,x) ' 

W>V = l + Y, m z_z iV f = ^(0 = ^AW™i(f) ■ (7) 

Here /ij = G^m^ where is the mass of the source. The geometric asymptotia is con- 
trolled by the completely homogeneous neutral energy density p n eutrai, which may involve 
a cosmological constant, usual neutral dark matter, baryons etc. The scale factor of the 
universe R(t) obeys the standard Friedmann equation H 2 = 8n ^ N p ncutral , with the charged 
dark matter contributions cancelled by the gauge repulsion terms thanks to q/m = ^/Gn 
for all the charged masses in the system. Indeed, this no-force condition between any pair 
of charged masses guarantees that the charged dark matter only contributes to local dis- 
tortions of the space-time geometry through fl(t,x), which is a harmonic function of the 
physical coordinates y = R(t)(x — Xi) . 

The solutions of the McVittie family in general FRW environs and their charged gen- 
eralizations have been questioned on the basis of the singular structures in the geometry 
[29, 30]. However, those issues only arise if one takes (7) to arbitrarily short distances. In 
the application we envision, the solutions (7) are naturally cut off by the finite size of the 
sources. As always, the cutoff scale is set by the inverse mass of the stable charged dark 
matter particle which sources (7), and which, by our discussion above, is bounded from 
below by £uv ~ ~ A -1 >, g~ l Mp\ [12]. Thus the solution (7) should only be trusted 
down to distances \y\ £ g~ x Mp\ , where the potential function Q is Q ~ 1 + (9(10~ 30 ), for 
the case of gravitationally coupled particles with m ~ TeV. At shorter distances, one would 
have to use the full description of the matter source, which would regulate the solution just 
as the interior solution regulates the field of a planet at distances shorter than the planet's 
radius [31]. In any case, this will occur far far away from any strong gravity regime near 
a singularity, where the gravitational fields are practically still in the linear regime, and so 
whatever the issues involving McVittie and singularities, they need not concern us here. 



Non-geodesic motion of charged particles 

Since different charged masses reside at fixed comoving coordinates Xi, their separation 
increases in time according to ^ = Hy, which is precisely the Hubble law with H set by 
only the neutral matter contributions. Indeed, if we consider at — const slice of (7) and 
evaluate the spatial distance between eg. the origin and a location of a charge at x — Xj, to 
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leading order we find it to be given by the proper distance between the origin and the outer 
surface of the particle, y = R\xi\ — £uv + A*i hi(i?|xi|/Ajv)> w here is the short-distance 
cutoff of the geometry near the charge. As time goes on, this distance increases according 
to ^ = H [y + 0(£\jy, Hi, Hi ln(y))) , with the time rate controlled by H of Eq. (7), 

FT - R FT 2 - 87VGN n (0\ 
ti — ti — Pneutral- [°) 

Note, that these trajectories are not geodesies: they involve the U(l)' Lorentz forces oc 
qi Fj ^ v which encode the gauge repulsion of the charge by all the other charges in the 
system. This is why they can stay put at X; L = const, and move relative to each other with 
recession velocities set up only by the neutral mass density. In other words, the charged 
mass density contributions cancel exactly in (7), because both the charged sources and the 
probes have the same q/m, which satisfies (q/m) 2 = Gn, as anticipated by Eq. (2). 

Geodesic motion of neutral particles 

On the other hand, neutral particles follow the geodesies of (7). The geodesic equations with 
the above notations read 



A2 



tt 2 



, F7 2 \7 '\ 2 (+'\ 2 

2 x') = ^^d £ Q 2 -^fd,Q- 2 , (9) 

(s^ = ^ ft n- 2 - ^ (i^n 2 ) , (io) 

+ R 2 Q 2 \x'\ 2 = -M 2 , (11) 



where the last equation follows from the afflne reparameterization invariance of the geodesies, 
and M. is the integration constant that gives the afnne parameter in the units of proper time. 
If we now follow the motion of a probe particle over sufficiently short times T < H -1 , such 
that R ~ const, then Eq. (10) implies that t' /Q 2 ~ const over the course of this observation. 
Picking this ratio to be t' /Q 2 ~ M. so that a particle infinitely far from the massive sources 
can remain at rest in comoving coordinates {i.e. ensuring that we get \x '\ — > as Vt — > 1), 
we obtain from Eq. (11) 

R 2 n 2 \x '\ 2 ~M 2 (n 2 -l)~2M 2 V. (12) 

to leading order in the harmonic sum V = J2i R\s-g y ^ e ra tio of R 2 Vl 2 \x '\ 2 to (t 1 ) 2 /VL 2 
gives the comoving kinetic energy per unit mass of a neutral particle, or equivalently, the 
square of the comoving velocity of a neutral particle, 



R 2 W\ 2 



~(l--L)~2V. (13) 



(t> 2 ) ~~ V ^2 

So the total velocity of a neutral particle is Wtotai = "^Hubble + V, where ^Hubble is the Hubble 

— * 

flow relative to our initial observer, and V is the total 'peculiar' velocity of the probe. 
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We now make some simplifying assumptions: (i) we ignore the asymmetries in the nearby 
mass distributions and treat the background as spherically symmetric about any observer's 
vantage point, and (ii) we view the probe's motion as being 'virialized' with all the sources 
acting on it, such that its motion corresponds to merely marginal binding with the masses 
surrounding the observer. The latter approximation randomizes the peculiar velocity of the 
probe, which is reasonable in order to extract the leading order bulk velocity. Indeed, the 
'average projection' of the net velocity on the direction of the Hubble flow relative to our 
initial observer, which due to our virialization assumption means that the observer gathered 
the information from many probes in a patch of the sky with random direction of their 
peculiar velocity, to leading order yields (^Hubble -V) = 0. Hence the net velocity squared of 
a particle at a distance R comparable to the horizon scale Hq 1 (i.e. that corresponds to the 
comoving distance |Ax| ~ 1) from an observer is 

l^totail 2 = I^HubbicI 2 + |V| 2 - H 2 R 2 + 2V . (14) 

Now imagine that the extra term in Eq. (14) is interpreted by the observer as a correction to 
the Hubble recession. It follows that a neutral particle, with q/m = 0, will 'see' an effective 
Hubble parameter H given by 

H o^H 2 + — . (15) 

Taking the continuum limit of the charged particle distribution and recalling that Hi = 
G^rrii, where is the mass of the individual charged dark matter particle, we see that 
^ = r3 2 |Ax| )| A ^| 1 — *■ ^y^Pchargcd, where as noted before we evaluate the potential at 
the comoving scale corresponding to the horizon distance, |Ax| ~ I. Substituting into (15) 
and using H 2 = Sn ^ N p ne utrai, as per Eq. (7), finally yields 

tti tti RkGn 8hGn 
~ 3 — Pchargcd = — 3~ Ptotal • ( ' 

This is precisely our starting Eq. (1) with (^) pro bc = 0, as it should be for the case of a 
neutral probe. 

For both charged and uncharged probes, we have now seen that the linearized limit of the 
exact solutions precisely confirms the Newtonian cosmology argument, yielding the expected 
contribution from gauge repulsion. One expects, on the basis of gauge symmetries of gravity 
and gauge theory, and the resulting energy and charge conservation, that the central features 
of this argument should remain valid beyond the linear limit (or more precisely, that there 
exists a gauge in which the theory effectively linearizes over the relevant scales, as for example 
in the derivation that led to (7)). So indeed the neutral particles feel the attraction of the 
total mass in the universe, while the charged ones do not sense the attraction of the charged 
component due to its repulsive gauge force. 

Local cosmology from the Newtonian gauge potential 

So far we have obtained our 'corrected' Friedmann equation (2) using intuitive arguments 
based on the motion of distant particles. Whereas the neutral particles follow geodesies, the 
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charged particles do not, since they suffer an extra acceleration due to the gauge repulsion. 
In the previous two subsections, we explained how this gives rise to 'probe dependence' in 
the Friedmann equation (2). To further allay any concerns that our 'corrected' Friedmann 
equation (2) has been derived 'too heuristically', we will now reinforce the derivation of Eq. 
(16) using an alternative approach, extracting it directly from the metric (7). 

Our aim now is to show that Eqs (15) and (16) follow from rewriting (7) in the 'New- 
tonian gauge' as in, for example, [36], and evaluating directly the Newtonian potentials at 
subhorizon distances. This is precisely consistent with our approach so far. So, let us begin 
with the full solution for the metric (7) and expand about an observer at the origin, so that 



ds 2 ~ - 



d&_ h 2f,A: | x i x .i 
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+ 0(x 3 ) , (17) 



where 



fig(t) =fi(t,0), ki(t) 



on 

dx i 



kij (t) 



d 2 n 



(t,0) 



dx l dxi 



(18) 



(t,o) 



Note that we do not need to worry about any of the gradients diverging since we have 
introduced a natural cut-off, £\jy, around the core of each particle. Clearly this metric may 
now be thought of as a local cosmological perturbation about the background Hubble flow 
set-up by the neutral matter content, p ne utrai(^)- 

In the absence of this perturbation, one can extract the background Hubble parameter 
by writing the FRW metric as a perturbation about the observer's local Lorentz frame on 
sub-horizon scales (see, for example, [36]). This is done by going to 'Newtonian gauge', 
where the Hubble parameter is given in terms of the Laplacian of the Newtonian potential, 
evaluated at the observer's position, H 2 = §V 2 ?/>|obs- Clearly, a perturbation about the 
cosmological background complicates things, but only slightly. To proceed with the same 
method of evaluating the Hubble parameter as in the absence of local inhomogeneities, we will 
again rewrite our perturbed solution (17) as a sub- horizon perturbation about the observer's 
local Lorentz frame, in Newtonian gauge. We begin with the coordinate transformation 
(t,x) — > (r,y), where 



t = T(t) 



x = 



■H(T)n,,(T)\y\ 



y 



l 



R(T)Qo(T) 
1 

+ R(Tyn (Ty 



l + l H 2 {T)\y\ 2 



1 



i?(T) 2 fi (T)3 
l,r, 
4 



(y-k(T))y--\y\ 2 k(T) 



(y- k(T)) 2 - -\k{T)\ 2 \y\ 2 ) y - Uy ■ k(T))\y\ 2 k(T) 



(19) 



+ C% 4 ) , (20) 



and T'{t) = Q. {T). This immediately puts the metric (17) in the desired form, 

ds 2 ~ -(1 + 2<f))dT 2 + (1 -2i))dy 2 , (21) 
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where the Newtonian potentials are given by 



<t> 



2R(T) 2 tt (Ty 



5(y • k{T)f - \k{T)?W 2 



y l y 3 h 3 {T) 

2R{T) 2 Q {Tf 



+ 0(y 3 ) , (22) 



i?(T) 2 fi (T) 4 



(y-k(T)) 2 -\\k(T)\ 2 \y\ 2 



2R(T) 2 tt (T) 3 



(23) 



These equations describe the gravitational field in the vicinity of the observer, as felt by 
neutral probes following geodesies of the metric (7). To extract the effective Hubble pa- 
rameter measured by the observer recording the recession of these probes, we simply take 
the Laplacian of the Newtonian potential, ip, as in [36], and evaluate it at the observer's 
location. This yields 



Hi 



^\y=6 = H 2 (T) 



+ 



(T,0) 



3R(T) 2 n (T) 3 3R(T) 2 n (T) 



(24) 



As expected, the background Hubble parameter gets shifted due to the charged dark matter 
components. To recover Eq. (16), we recall that Sl = l + V, where V = Ei ms-i^ and 
rrii is the mass of the dark matter particle at Xj. Taking the continuum limit, the potential 
V becomes an integral V — > / d 3 x' ^j^zfrj , where p(x') is the local comoving density of the 
charged dark matter. It follows immediately that 



- d A x 



, 3 ,G N p(x')x' 



R(T)\x'\ 



(T,0) 



-47TG 



N 



p(Q) 

R(T) 



(25) 



We now make use of our simplifying assumptions from before, and ignore asymmetries in 
the nearby mass distributions, treating them as spherically symmetric about the observer. 
This eliminates the first derivative terms, V x fiL T g\ — > (which appear non-linearly in Eq. 
(24) anyway, and hence really correspond to higher multipole corrections in the expansion, 
and are essentially negligible for our purposes), so that 



Hi 



H 2 (T) + 



8ttG 



N 



( P®_ 



3 \R(T) 3 tt (Ty 



H\r) + 



8ttG 



N 



Pcharged(r) (1 + 0(V)) , (26) 



where p c har g ed(T) = p(0)/-R(t) 3 is the physical energy density of the charged dark matter, 
to leading order. As stated earlier, the 0(V) corrections are never larger than 10~ 30 . Thus 
we recover Eq. (16), providing a rigorous check of our modified Friedmann equation (1) for 
neutral probes. 
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Validity of the multi- centered black hole solution 

So far, we have treated the distribution of the charged dark matter particles as essentially 
static in comoving coordinates, making use of the multi-centered black hole solution (7). 
That is of course unrealistic: the dark matter particles won't just stay put at some fixed 
values of Xi, but will move around with peculiar velocities, as large as (9(100) km/s. Due 
to the different boost properties of gauge and gravitational conserved currents (the latter 
being of course the conserved stress energy tensor), the no force condition between different 
charged particles will be violated, leading to net forces between them that are proportional 
to their relative velocity. Thus the whole array (7) will really be dynamical, being unstable 
to nonzero velocity perturbations, and with individual constituents forever bouncing around 
in each other's fields [32, 33]. However, this needn't concern us too much here, as we will 
now explain. 

First of all, the dark matter particles are packed together rather closely: since their 
number density in the universe is n ~ p C r/ m ~ 10~ 12 eV 4 /m, their interstitial distance is 
£ DM ~ 1/n 1 / 3 ~ (m/eV) 1 / 3 mm. If the dominant component of dark matter is a TeV- 
mass WIMP, there is a dark matter particle in every cube of side ten meters or so in the 
extragalactic medium. The density of charged dark matter in extragalactic space could be 
more dilute. However if the mass of charged dark matter is also in the WIMP range, which 
would follow from the conjecture of [12] and our taking the U (1)' strength to be gravitational, 
their separation will also be within an order of magnitude of WIMP's, i.e. ~ (9(10) meters. 
On the other hand, if they are black holes, with masses m > M P \ ~ 10 18 GeV, their interstitial 
distance in the extragalactic medium will be of the order 1000 kilometers. In any case, in the 
units of the Hubble scale, they are very close together. Clearly, their separation would be 
yet smaller in condensed structures, such as galaxies, where the matter density may exceed 
the uniform extragalactic density by a billion times and more, being p ga i a xy ~ 10~ 3 eV 4 or so. 
Thus, when considering the dark matter interactions, to leading order we can safely treat 
them as swarms of particles in flat space. 

In this case, the analysis of [33] shows that they will affect each others trajectories, owing 
to the fact that gravity has infinite range. However, they will gravitationally coalesce only 
when the impact parameters are of the order of their horizons. If they are WIMP-scale stable 
particles, the dynamics of the theory which completes their description in the UV will become 
important at much longer scales, £uv ~ m ~ 1 as P er the conjecture of [12]. Even on the rare 
occasions where they form black holes with masses M ^> m, the weak gravity conjecture 
of [12] suggests that such black holes will decay back into the charged constituents, being 
merely metastable transients. Generically, however, such dark matter particles will simply 
zip around each other, mostly missing each other. In any case, at very large distances, to 
observers moving away due to Hubble recession, they will appear as nearly static swarms of 
mass and charge, whose conservation will justify using (7) as the appropriate leading order 
approximation 7 . Indeed, as long as the distance to the observer exceeds L ~ v peculiar /H, 

7 In the papers [34] the statistical mechanical description of the gas of charged black holes has been studied, 
with further interesting possibilities for cosmological applications found to occur if the gas condenses. Our 
analysis as exhibited in this work is appropriate for the dilute gas limit. Here, we will not address the short 
interstitial distance approximation where condensation may occur because of possible model dependence of 
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such that the Hubble velocity is the dominant driver of the change in relative separation, it 
is perfectly reasonable to make use of the multi-centered black hole solution (7). 



4 Discussion 

We have clearly shown that the presence of dark matter charged under a dark U(l)' can 
lead to a mismatch in the estimates of the dark matter density, depending on whether we 
use neutral or charged probes. Observers who fail to take into account the gauge repulsion 
between the galaxies which house type la SNe, and may contain charged dark matter, will 
misinterpret the enhanced cosmic recession as an otherwise mysterious 'anti gravity' effect, 
inferring a smaller equation of state for dark energy than is actually the case. For example, 
if dark energy really were a cosmological constant, they may erroneously deduce that the 
dark energy equation of state would lie beyond the phantom divide, w < — 1. Although we 
have focussed on the role of probes used in observations, our corrected Friedmann equation 
(2) includes the possibility that the background charge-to-mass ratio can also be redshift 
dependent. This can occur because more charged dark matter will have been attracted to 
collapsed structures at later times. 

In the presence of net charge and associated gravitational strength long range forces there 
can be other effects. Let us comment on one very interesting effect that may also arise from 
inhomogeneities in our model: local variations in the density of charged dark matter can 
lead to local variations in the Hubble flow. This could be relevant for explaining the recent 
claims of [17]. Suppose that in some region of space, of size considerably smaller than the 
Hubble length, but big enough to accommodate many galaxies, there is an overdensity of 
charged particles bp. The presence of extra charged matter will lead to an extra repulsive 
force acting on the particles in that region, as compared to the forces acting on neutral 
particles, directed radially outward from the center of the overdensity. The velocity of the 
charged probes inside this region will therefore be larger than the mean recession velocity 
estimated by the Friedmann equation (2). To get the scale of the effect, we can proceed 
as follows. From Eq. (13) we see that the charged probes will move a bit faster than the 
neutral ones, by an amount 

Av ~ ^ VW . (27) 

Now compare this to the mean recession velocity in a region of space of size R, given by 
v ~ HR. Assuming that inside the overdensity the charged particles are distributed ap- 
proximately uniformly, with density <5p c har g ed = £ Per, where p cr is the critical density, we 
find 

Av2 = ^! ^ 6V „ Charged = , (2g) 

v 2 H 2 8irG N p cr R 2 p cr 

This means that the excess velocity will scale as the square root of the galactic charged mass 
density, in units of the critical density. 



the cutoff effects, that may come in because of the weak gravity conjecture of [12], on the results of [34]. 
This remains an interesting question to be considered elsewhere. 
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In the recent paper [17] it has been reported that there may be evidence for anomalous 
flows in an observed region of the universe of a typical scale ~ (9(100) MPc, out to about 
300 MPc. Since v 2 ~ ^-G N p cr R 2 ~ H 2 R 2 , its numerical value at the scale of R ~ 100 - 300 
MPc, which is about 0.01 to 0.03 of the Hubble length, ranges between v 2 ~ 10~ 4 and 10~ 3 . 
Thus, the relevant excess velocity is given by 



On the other hand, [17] report excess velocities of the order of 1000 km/sec, which, in units 
where c = 1, yields 0.003. This gives a value for Av 2 of about 10~ 5 . To reproduce this scale, 
the charge overdensity £ need only be in the range of 10 percent to a percent, diminishing 
as we move away from its center. This appears perfectly plausible, allowing for a local 
interpretation of the claims of [17], if they are proven to be a real physical effect, and 
without any need of invoking superhorizon physics. It may therefore warrant a closer study 
and comparison with the data. 

Further, one may also expect that to subleading order, the CMB peaks could be affected 
because the charges will alter the formation of structure. The charges of the same sign will 
lead to a suppression os structure formation, since they will 'screen' part of the collapsed 
structure's mass and impede its further growth. This can shift the CMB peaks some. To 
get an idea by how much, we note that if £ is the fraction of charged dark matter, the 
correction in the suppression of structure must go as £ 2 , just because the extra force requires 
two charges to operate, if we taking the initial density perturbation is adiabatic. In this case, 
the suppression will only occur as the additional charged dark matter particle are pushed 
away by the charges already inside the structure. This is also evident in our equations where 
the effects of charge go as (q/m) 2 . Hence, small amounts of charge will not affect CMB a 
lot. 

One should also worry about the possible effects of the extra U(l)' on nucleosynthesis. 
If too much energy in the early universe is deposited in the massless U(l)' mediators, it 
could compete with the photon energy during nucleosythesis, spoiling light element abun- 
dances. Indeed, we know that the extra light matter contributions to the expansion rate of 
the universe during nucleosythesis must be below photons by at least about two orders of 
magnitude. The question then is, how much U(l)' will be produced early on? Clearly the 
answer is model-dependent; however, generically due to the weakness of the U(l)' couplings, 
one expects the production of U(l)' to be suppressed. Suppose, that the net charge in the 
universe is produced by the decay of an inflaton condensate that was charged under U(l)'. 
There will be no significant U(l)' particle production during inflation itself, due to the con- 
formal invariance of any U(l) sector, well known from the example of Maxwell's theory in 
inflation. At the end of inflation, however, when the inflaton starts to oscillate around its 
minimum, there may be particle production due to the mflaton-£7(l)' coupling ~ gA^J^, 
where g is the U(l)' coupling and is the inflaton current. If the spatial fluctuations of 
the inflaton produce a net U(l)' dipole over a Hubble volume at the end of inflation, the 
total U(l)' power will be of the order of P ~ o; 4 |p| 2 ~ H 4 x q 2 /H 2 , since the characteristic 
frequency is given by the inflaton mass, oo ~ m ~ H, and the separation of charge in the 
dipole moment is maximized by the Hubble scale, / ~ 1 / iif , so that p ~ q/H. The energy 
density deposited in U(l)' mediators will be pu(iy ~ PTh/V h ~ q 2 H 4 . The net charge in 




(29) 
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the Hubble volume will be comparable to its net mass M# ~ pVn ~ (GV-ff) -1 , and so q <, 1, 
implying that pu(iy ~ -f^ 4 <C Preheating, because it is comparable to the gravitational particle 
production which is negligible when inflation occurs at, or below, the GUT scale. Similarly, 
the repopulation of the universe by the U(l)' mediators, that can occur by the collisional 
emissions of the U(l)' quanta, will be suppressed relative to the photons by the square of 
the ratio of their couplings, ~ g 2 , and it will be negligible for gravitational couplings. 

Thus, the most interesting question one can ask here is, just how much charge could there 
be? The answer should involve a detailed survey of various cosmological phenomena, and to 
have it we need to know the precise distribution of charges in the universe. Evidently, it is 
impossible to write down the formula for (q/m)(z) without a detailed study of how structures 
form, and acquire charged dark matter. In the presence of charged dark matter, however, 
it must be non-zero for low red-shift observations like type la SNe, simply because when 
structures form by gravitational instability, some of the charged dark matter will fall in along 
with neutral matter. Thus, at least in the linear regime, one expects q/m to be comparable 
to the total charged dark matter fraction at low redshift. Before this epoch, at large redshift, 
the charged particles will be distributed approximately uniformly in the universe. If during 
that epoch we use neutral probes, like CMB as observed by WMAP, which involve only 
the propagation of photons of zero charge-per-mass in the background geometry set up by 
total masses, to the leading order they will not directly discern the presence of the charged 
particles. Clearly, at subleading order there will be effects, since the charges affect mass 
distributions which in turn affect the CMB anisotropics. Of course, to be able to estimate 
these effects, one should i) produce a precise parameterization of (q e ff/m)(z) for a given 
class of cosmological probes, beyond the simple 'step function' distribution sketched herein, 
and ii) account for possible systematic effects involving various uncertainties in present data. 
This task is beyond the scope of the present work, and we hope to return to it in the future. 
Here, our aim is to merely indicate that the search for anomalies in the velocity fields of the 
universe, cosmic structures and dark energy equation of state might reveal the presence of a 
new long range force, or at least place stronger bounds on them. 

To this end, we can still get plausible numerical estimates of the scale of the effect. 
Imagine for simplicity that dark energy is a cosmological constant, with equation of state 
wde real = — 1- Imagine also that the density of charged dark matter is, say, 5% of the real 
dark energy density, and take the U(l)' strength to be exactly gravitational, assuming that 
the galaxies accreted charged dark matter in the same proportion to neutral matter, as in 
extragalactic space. From Eq. (5) one then gets that u>de eff = —1-053. This is the dark 
energy equation of state inferred by an observer who does not know of the presence of dark 
charges, and correlates the low redshift supernova data as if all the large scale recession 
were controlled by gravity alone. As it stands, this might not appear as much. Yet note 
that it is near the limits of observationally testable deviations of dark energy equation of 
state from w — — 1, and that it demonstrates how the data may suggest the crossing of 
the phantom divide without any exotic new physics such as ghosts or similar pathologies. 
Actually, the effect may be even larger if the dark U(l)' is in fact stronger than gravity, such 
that (q/m) 2 > G^. This is consistent with the 'gravity as the weakest force' conjecture of 
[12], and may even be favored by it. In such a case, clearly, one would expect to have smaller 
^charged in order to avoid disturbing too much the standard picture of structure formation, 
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however this may be compensated by the increase in Gj^(q/m) 2 . To get a more accurate 
idea of the possible range of wde eff, one would have to undertake a detailed study of bounds 
on charged dark matter from structure formation and other cosmological observations. 

Note, that in a more accurate accounting for the effects of the charged dark matter 
population, one should also consider the influence of inhomogeneities in the cosmological 
backgrounds on the propagation of light. The inhomogeneities alone are known to affect the 
luminosity-distance relationship, perturbing the expression for angular diameter distance as 
a function of z. In a careful analysis this ought to be properly included [35]. In the presence 
of charged dark matter, the inhomogeneities can arise from both neutral and charged particle 
density fluctuations, and should also be carefully considered. We will ignore the detailed 
derivation of those effects here, keeping in mind that it can be done such that its influence 
is subtracted out of (3) when considering the physical effects of dark charges. We hope to 
return to this issue more carefully in the future. 

In closing, let us summarize our findings. If there are dark matter particles which are 
charged under an unbroken U(l)', of gravitational strength, they may have very interesting 
long range effects. Such states may arise in string theory, as the heavy BPS states. In the 
exact supersymmetric limit, their charges are equal to their mass, and, for states with masses 
close to the supersymmetry breaking scale, remain comparable even after the breaking 8 . We 
stress, that for the phenomena described here, the only required condition is that the charges, 
and the forces they source, are comparable to the particles' masses and their gravitational 
attraction, respectively; they do not need to be the same. Now, when their masses are of 
the order of string scale or greater, they are interpreted as extremal black holes. However, 
it has been conjectured [12] that in any UV complete framework gravity should be the 
weakest of forces, which implies a low cutoff in the gauge theory and the existence of stable 
states lighter than this cutoff. In any case, the long range fields of the charged dark matter 
can be described by the extremal black hole geometries below the cutoff, or equivalently 
outside of the sources. For particles with the same sign charges, their gauge repulsion will 
compensate gravitational attraction, weakening their net binding. If they partake in the 
matter population of the universe, measuring the masses with orbits of different probes, as 
is typically done in astronomy, may yield a mismatch between effective masses of cosmic 
objects. Specifically, we have seen that the flows of dark matter structures, such as galaxies, 
may yield smaller values if the gauge repulsion is not accounted for. In particular, the 
observation of type la SNe which reside in galaxies that contain charged dark particles, may 
misinterpret their faster cosmic recession due to weaker forces for an extra 'antigravity', 
simulating effective dark energy equation of state smaller than the real one, and in some 
cases mimicking w < — 1. We have also noted that in the presence of charge overdensities, 
there may be local variations of intergalactic forces, that can give a larger 'Hubble Flow' 
in those regions of space. Such phenomena can be confused with a dynamical dark energy, 
or superhorizon effects. Thus it is important that such effects are taken into account when 
attempting to extract the properties of dark energy from cosmic motion and geometry, as 
they may be contaminated by the long range forces in the dark sector. 

Corrections to the coupling will be negligible in the extreme weak coupling when U(l)' is unbroken. 
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